Single-spin azimuthal asymmetry in exclusive electro-production of 7r + mesons on 

transversely polarized protons 
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Exclusive electroproduction of tt + mesons was studied by scattering 27.6 GeV positrons or elec- 
trons off a transversely polarized hydrogen target. The single-spin azimuthal asymmetry with respect 
to target polarization was measured as a function of the Mandelstam variable t, the Bjorken scaling 
variable xb, and the virtuality Q 2 of the exchanged photon. The extracted Fourier components of 
the asymmetry were found to be consistent with zero, except one that was found to be large and 
that involves interference of contributions from longitudinal and transverse virtual photons. 
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PACS numbers: 13.60.-r, 13.60.Le, 13.85.Lg, 14.20.Dh, 14.40.Aq 



Generalized Parton Distributions (GPDs) [H, [2, 3] pro- 
vide a three-dimensional representation of the nucleon 
structure at the partonic level correlating the longitu- 
dinal momentum fraction of a_parton with its trans- 
verse spatial coordinates 0, [a S 0] ■ The possibility 
to study GPDs relies on factorization theorems proven 
in the framework of perturbative quantum chromody- 
namics for hard exclusive processes at leading twist, in 
particular for hard production of mesons by longitudi- 
nal virtual photons Q . For recent theoretical reviews, 
see [13, EL [ll | ■ 

In the description of hard exclusive electroproduc- 
tion of pseudoscalar mesons at leading twist, only the 
two GPDs H and E appear. Spin-averaged and spin- 
dependent cross sections are sensitive to different com- 
binations of H and E. It was predicted that for ex- 
clusive production of 7r + mesons on transversely polar- 
ized protons by longitudinal virtual photons the interfer- 
ence between the pseudovector (oc H) and pseudoscalar 
(oc E) contributions to the cross section leads to a large 
proton-spin related azimuthal asymmetry [I! El. Un- 
like the spin-averaged cross section, this asymmetry is 
directly proportional to the sine of the relative phase be- 
tween H and E. It was shown that next-to-leading or- 
der corrections in the strong-coupling constant a s cancel 
in the asymmetry [HI, [3]. No GPD-based model pre- 
dictions are available for the production of tt + mesons 
by transverse virtual photons as no factorization theo- 
rems exist for this case, but also because the leading- 
twist contribution is expected to be dominant. Measure- 
ments of the asymmetry are considered to be a valuable 
source of information about possible contributions from 
transverse virtual photons [17] . In a Fourier expansion 
of the proton-spin-dependent part of the hard exclusive 
pion electroproduction cross section 18( the only leading- 
twist contribution to the asymmetry from longitudinal 
virtual photons is the sin(<^ — 4>g) Fourier amplitude, 
which can be used to test GPD models. All other ampli- 
tudes involve contributions from transverse virtual pho- 
tons. Here, following the Trento conventions [l9j . <j) and 
4>s are the azimuthal angles in the proton rest frame of 
the pion-momentum and the proton-polarization vectors, 
respectively, measured about the virtual-photon momen- 
tum vector relative to the lepton scattering plane. For 
recent theoretical analyses of exclusive pion electropro- 
duction, see 17, 2(| 21 1 . 

The Hermes collaboration has previously performed 
measurements of the spin-averaged cross section 22| and 
the single-spin azimuthal asymmetry in exclusive 7r + elec- 
troproduction on longitudinally polarized protons [23l |. 
This letter reports the first measurement of the single- 
spin azimuthal asymmetry for the hard exclusive reac- 
tion ep^ — > enir + on transversely polarized protons. The 



kinematic variables relevant for the analysis of this pro- 
cess are the squared four-momentum of the exchanged 
virtual photon q 2 = —Q 2 , the Bjorken variable xb = 
Q 2 / (2M p v), and the squared four-momentum transfer 



(q - Pir+) 2 - 



Here, M p is the proton mass, v the 



t 

energy of the virtual photon in the target rest frame, 
and Pk+ the four-momentum of the pion. Instead of i, 
the quantity t' = t — to is used in the analysis, where —to 
represents the minimum value of — t for a given value of 
Q 2 and xb- 

The data corresponding to an integrated luminos- 
ity of 0.2 fb -1 were collected with the Hermes spec- 
trometer in the years 2002-2005. The 27.6 GeV 
positron or electron beam was scattered off the trans- 
versely nuclear-polarized gaseous hydrogen target inter- 
nal to the Hera storage ring at Desy. The open- 
ended target cell was fed by an atomic-beam source [HI 
based on Stern-Gcrlach separation combined with ra- 
diofrequency transitions of hydrogen hyperfine states. 
The nuclear polarization of the atoms was flipped at 
1-3 minute time intervals, while both this polarization 
and the atomic fraction inside the target cell were con- 
tinuously measured [26| . The average magnitude of the 
transverse polarization of the target with respect to the 
beam direction was |Pt| = 0.72 ± 0.06. 

Events were selected with exactly two tracks of charged 
particles: a lepton and a pion. Furthermore, it was 
required that no additional energy deposition was de- 
tected in the electromagnetic calorimeter. The Her- 
mes geometrical acceptance of ±170mrad horizontally 
and ±(40-140) mrad vertically resulted in detected scat- 
tering angles ranging from 40 mrad to 220 mrad. Lep- 
tons were identified with an average efficiency of 98% 
and a hadron contamination of less than 1% by using an 
electromagnetic calorimeter, a transition-radiation detec- 
tor, a preshower scintillation counter, and a dual-radiator 
ring imaging Cerenkov detector [27l |. Pions were identi- 
fied in the momentum range 2 GeV < p < 15 GeV using 
the Cerenkov detector. For this momentum range the 
pion identification efficiency was on average 99% and the 
contamination from other hadrons less than 2%. The 
kinematic requirement Q 2 > 1 GeV 2 was imposed on the 
scattered lepton in order to select the hard scattering 
regime. 

The single-spin asymmetry for exclusive 7r + produc- 
tion with unpolarized (U) beam and target polarization 
transverse (T) to the lepton {€} beam direction is defined 
as 

a I a a \ 1 do " T (<Ms) ~ do-lp, 4> s ) m 
AvrAMs) - |^ daWs) + daWsr W 

where dcr T(x) (0, cf>s) = dcr U u(</ ) ) + Pr d<ruT,*0A> fa) is a 
sum of the spin-averaged and spin-dependent cross sec- 
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tions, with Pt/|Pt| equal to 1 (—1) for the f (J.) orien- 
tations of the transverse target polarization vector Pt\ 
Both numerator and denominator of (Q]) can be Fourier- 
decomposed [lH, respectively, as 



dcruT/(0, 4>s) oc 2(sin(</> - <l>s)}vT,e sin(0 - <^>s) + 



(2) 



where the ellipsis denotes five more terms omitted here 
for brevity, and 



dtruu(^) oc 1 + 2(cos^)uu cos<fi 

+ 2(cos(20))uu cos(2<^). 



(3) 



Ideally, the Fourier amplitudes in ([2]) , which provide most 
direct access to the photoabsorption subprocesses, should 
be measured, e.g., 

(sm((f) - 0s))uT,£ 

__ / <j<j>dc/)s sin(0 - <ps) dtr UT ^(0, (j>s) ( 4 ) 
Jd4>d4> s dcr uu (0) 

For experimental reasons, mainly to minimize effects of 
the Hermes spectrometer acceptance in </>, the Fourier 
amplitudes associated with the asymmetry ([l} were ex- 
tracted instead, e.g., 



1 UT,i 



>>-<t>s) 



4tt 2 



s sm(</> - <ps) 



do-uu(^) 



(5) 



Similar equations hold for the other five amplitudes. 
These amplitudes embody all the essential information 
that could also be extracted from ([2]). For small (or zero) 
values of (cos </>)uu and (cos(20))uu, the amplitude in (O 
corresponds to the one in ((4]). 

The set of six Fourier amplitudes of the asymmetry 
was obtained from the observed 7r + event sample using a 
maximum likelihood fit. The distribution of events was 
parameterized by the probability density function Af pal 
defined as 

■A/" P ar(Pr,^,^s;»7uT/) = 1 + Pt*4utv(</>, <As;»7ut/)> 

(6) 



where 



-4uT,f(</>, 0s;»7ut,£) 

= A^*~ M sin(^ - S ) + 4S +0S) sin(0 + <f> s ) 

+ sin <p s + 4t (2 /" 0s) sin(20 - cf> s ) 

+ A^f~* a) sm(3<£ - <fe) + Af T ( f +4,s) sin(20 + <f> s ). 

(7) 

Here, t]tj T e represents the set of six Fourier amplitudes 
of the sine-modulation terms in (T5|). 



Within the maximum likelihood scheme [28| , the loga- 
rithm of the likelihood function to be minimized is taken 



as 



N_ 



^ln[l + P^ UT , 



(8) 



i=l 



where N-+ = N [ + 



N l + is the total number of events 



in the selected data sample, and 

The raw results from the likelihood minimization of ([5]) 
were corrected for background contributions in order to 
estimate the true results for exclusive 7r + production: 



-4,, 



A T - bA h 
1-6 



(9) 



Here, A T stands for one of the six Fourier amplitudes in 
^ut l ( see ©i ©): b and A\> for the fractional contri- 
bution and corresponding Fourier amplitude of the back- 
ground, and A t for the resulting true amplitude. The 
background fraction is 



7T + 



cxcl 



(10) 



where N c + is the number of exclusive events in the se- 
lected data sample. 

The following analysis was performed to estimate the 
quantities in |9]). As the recoiling neutron in the process 
ep^ — ► enn + was not detected, the sample of "exclusive" 
events was selected by requiring that the squared missing 
mass M\ of the reaction ep^ — > eir + X corresponds to 
the squared neutron mass M% . The exclusive ir + channel 
could not be completely separated from the channels with 
final states ir + + X (defined as background channels for 
X =/= n) in which the tt + originates, e.g., from neutral- 
meson (mainly p°) decays, semi-inclusive processes, or 
nucleon resonance production, as their M\ values were 
smeared into the region around M 2 due to the experimen- 
tal resolution. These background events were subtracted 
from N v + following the method briefly outlined below, 
and previously employed in the analysis of the exclusive 
7r + cross section [221 ] . The exclusive ir + yield was ob- 
tained by subtracting the yield difference (N 7I + — N v - ) 
of the Pythia Monte Carlo simulation from that 
of the data, with both differences being independently 
absolutely normalized: 



N 



cxcl 



\data 



\ Pythia 



(11) 



The Pythia generator was used in conjunction with a 
set of Jetset [3(1 fragmentation parameters that had 
previously been adjusted to reproduce exclusive vector 
meson production data [3l| and multiplicity distribu- 



tions 



32| observed by Hermes. Exclusive production 
of single pions is absent in Pythia. Note that exclu- 
sive 7T _ mesons cannot be produced on protons. The 
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constraint on the invariant mass of the initial photon- 
nucleon system W 2 > 10 GeV 2 was applied, and the pion 
momentum was required to satisfy 7 GeV < p < 15 GeV. 
Both conditions, applied to the data and the Pythia 
yields, allowed for a better description of the data by the 
Pythia Monte Carlo simulation for values of M x out- 
side the region corresponding to exclusive n + production. 
The resulting M x distribution of iV™ 1 and its resolution 
of 0.7 GeV 2 were found to be consistent with that of a 
Monte Carlo sample of exclusive 7r + events normalized 
to the data (including radiative effects) [2^ |. 

An "exclusive region" in M x was defined by requiring 
—1.2 GeV 2 < M| < 1.2 GeV 2 . The lower limit corre- 
sponds to three times the resolution of M x , while the 
upper limit was set in order to minimize the (quadrati- 
cally) combined statistical and systematic uncertainties 
of the extracted Fourier amplitudes. A relative system- 
atic uncertainty of 20% was assigned to 7V e + cl , which cor- 
responds to the largest data-to-PYTHiA discrepancy out- 
side of the exclusive region [22] ]. As the M x spectrum 
of the positron-beam data is found to be shifted by ap- 
proximately 0.16 GeV 2 towards higher values relative to 
that of the electron-beam data, the exclusive region for 
the positron data is shifted accordingly. One quarter of 
the effect of this shift on the results presented below is 
assigned as a contribution to the systematic uncertainty. 

The values of A r and b in (|9|) are measured in the exclu- 
sive region. As the background originates from resolution 
smearing of events occurring at higher missing mass, 
in |9]) was assumed to be equal to the Fourier ampli- 
tude measured in the M x region between 1.9 GeV 2 and 
3.3 GeV 2 where the contribution of exclusive 7r + events is 
negligible. In that region A\> was found to vary smoothly, 
with values smaller than ±0.1, except for the sin^g mod- 
ulation for which it amounts on average to (0.25 ± 0.04). 
In order to account for a possible variation of A^ with 
M\ in the exclusive region, one half of the difference 
between A t and A r is conservatively assigned as a con- 
tribution to the systematic uncertainty of A t . 

The values of t 1 were calculated from the measure- 
ment of the four-momenta of the scattered lepton and 
produced pion by setting Mx — M n , which improved 
the i'-resolution by a factor of two. The kinematic 
range that contains the events used in the subsequent 
analysis is defined by the following requirements on the 
variables: -t' < 0.7 GeV 2 , 0.03 < x B < 0.35, and 
lGeV 2 < Q 2 < 10 GeV 2 . The mean W 2 value of the 
data is 16 GeV 2 . 

The dominant sources of systematic uncertainty are 
associated with the background subtraction and correc- 
tion, and the observed relative shift of the M x distribu- 
tions between positron and electron data. The contribu- 
tions due to the residual beam polarization of 0.02±0.03, 
the corresponding beam-spin asymmetry [23| . and the 
charged-track curvature in the transverse field of the 
target magnet, are found to be negligible. All these 



contributions, except for the target polarization scale 
uncertainty of 8.2%, are added in quadrature to yield 
the total systematic uncertainty. In addition, an esti- 
mate of the combined contribution to the experimental 
uncertainty from resolution smearing, acceptance, kine- 
matic bin width, and effects from the detector alignment 
with respect to the beam is determined u sing Monte 
Carlo simulation based on the GPD model [17j for the 
sin(</> — <j>s) Fourier amplitude only. The difference be- 
tween the amplitude extracted from the Monte Carlo 
sample and the corresponding model prediction calcu- 
lated at the average kinematic values of the Monte Carlo 
sample is added in quadrature to the total systematic 
uncertainty of ^tyxt ^ S ^- The largest experimental un- 
certainties are those due to detector acceptance and kine- 
matic bin width, and the determination of the target po- 
larization. 

Figure Q] shows the extracted Fourier amplitudes as 
a function of — t', xb, and Q 2 . For this measurement 
the average values of the kinematic variables are (— t') = 
0.18 GeV 2 , (x B ) = 0.13, and (Q 2 )=2.38 GeV 2 . The back- 
ground fraction b varies between (54 ±6)% and (62 ±5)% 
in the various kinematic bins. As xb and (Q 2 ) are corre- 
lated the average values of Q 2 vary in the four xb bins, 
namely, (Q 2 ) = 1.24, 1.57, 2.24, 3.91 GeV 2 . Analogously, 
the average values of xb vary in the four Q 2 bins, (xb) — 
0.07, 0.11, 0.15, 0.23. A separation of the contributions 
from longitudinal and transverse virtual photons to the 
Fourier amplitudes was not possible without measure- 
ments with different beam energies. 

The six Fourier amplitudes shown in Fig. [T] correspond 
to the following combinations of photoabsorption cross 
sections and interference terms cr^ n for photon helicities 
m,n = 0,±1 and proton-spin projections i,j = ±5 [Hj]: 

(sin(0 - <fis))vT,e oc [cosfllm (cr|+ + sa^ ~) 

+ - sin <V e(l + (0+0" - o+o~)]> 

(12) 

(sin(0 + <j>s))uT,t « b cos6>£lmer+I 

+ - sin <Ve(l + £)Im(CT++ - )], 

(13) 



(sin</>s) UT ,f oc [cos 9 ^Js{l+~e) Ima| ], (14) 



(sin(20- <t>s))wT,t « [cos6y/e{\ + e)Imcr +c [ 

1 (15) 
+— sin 9 elm o+_] , 

(sin(30 - (f>s))vT,£ oc cosG e1ma+±], (16) 
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FIG. 1: The set of six Fourier amplitudes (Aut,<) describing 
the sine modulations of the single-spin azimuthal asymmetry 
for unpolarized (U) beam and transverse (T) target polariza- 
tion, for the exclusive event sample. The error bars (bands) 
represent the statistical (systematic) uncertainties. The re- 
sults receive an additional 8.2% scale uncertainty correspond- 
ing to the target polarization uncertainty. 



(sin(20 + 0s)) UT ,£« [-sin0£lm<7+l], (17) 

where e is the virtual-photon polarization parameter, 
and 9 is the angle between the beam and the virtual- 
photon direction. Note that in the analysis presented 
here there is an integration over a range in 8, with 
cos# w 1 and 0.04 < sin# < 0.15. At leading twist, only 
(sin(</> — cj)s))\JT,e receives a contribution from only longi- 
tudinal virtual photons via <Jq^~ , while the other Fourier 
amplitudes are expected to be suppressed [§| by at least 
one power of 1/Q due to interference between contribu- 
tions from longitudinal and transverse virtual photons, 
and by 1/Q 2 due to terms involving only transverse vir- 
tual photons. 

Most of the Fourier amplitudes shown in Fig. [1] are 



small or consistent with zero, except A^ff . This am- 
plitude is found to be large and positive indicating a sig- 
nificant contribution from the transverse-to-longitudinal 
helicity transition of the virtual photon, i.e., 



(18) 



M* a+++ .Mo+o- + Mq_ ++ Mq-o- 



where .M M v M „ are helicity amplitudes with \j! [p) and 
v' (y) denoting the helicities of the pion (virtual pho- 
ton) and the neutron (proton) , respectively. These ampli- 
tudes are proportional to V~ & " +I/ ' . In the frame- 
work of GPDs, the amplitude Aio ^+ is associated at 

leading twist with virtual-photon helicity flip in the t- 
channel [la ], which is proportional to \f—t' and hence is 
expected to vanish for —t' — > 0. However, among higher- 
twist contributions the one that involves the parton- 
hclicity-flip GPDs Ht and Ht need not vanish at small 
values of \t'\. Moreover, in the more general framework 
of helicity amplitudes and the Regge model, ^4 UT ^ S re ~ 
ceives contributions from natural and unnatural-parity 
exchange [It], [33| , which allow it to remain constant as a 
function of — t', as the data in Fig. Q] suggest. Lack of pa- 
rameterizations of cr^ in involving transverse virtual pho- 
tons does not allow further interpretation of the corre- 
sponding Fourier amplitudes. Any model that describes 
exclusive pion production will need to describe not only 
the leading-twist Fourier amplitude, but also the other 
contributions to the target-spin azimuthal asymmetry. 
Of special interest in the present measurement is the 



Fourier amplitude A 



sin(0— 0s) 
UTV 



in case of production by 



longitudinal photons, which can be compared with GPD 
models. It is related to the parton-helicity-conserving 
part of the scattering process and is sensitive to the in- 
terference between Ti and £ HI, lij |: 



A 



UT,, 



M„ 



(i - ew 



4Mj ° 



2£ 2 Re(£*H) 



(19) 



where the transition form factors H. and £ denote con- 
volutions of hard scattering kernels and the pion distri- 
bution amplitude with the GPDs H and E, respectively. 
Note that in the models described below terms propor- 
tional to the coscj) and cos(20) modulation of the spin- 
averaged cross section are not included. In the measure- 
ment presented here these terms are not known, although 
they nonetheless contribute to the values of the extracted 
Fourier amplitudes. 

Figure [2] shows the extracted Fourier amplitude 



A 



UT,. 



as a function of —t' in more detail. The solid 
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and dotted curves represent the leading- twist, leading- 
order in a s calculations of this amplitude for longitudi- 
nal virtual photons using two variants of the GPD model 
of 20]. The modelling of the GPD E relies here, even 
at larger values of — i, on the dominance of the pion pole 
l/(mj — t) in the pion exchange amplitude, with the 
pion mass. Then £ is real and positive, and the value of 



.4' 



sin(0-0 s ) 
UT\. 



is typically predicted to be large and negative, 
while it must sharply vanish at the kinematic boundary 
— t' = (see solid curve). The data qualitatively disagree 
with such a simplified GPD model. The "Regge-ized" 
variant of the GPD-E model [20| . containing more than 
only a pion i-channel exchange, results in the dash-dotted 
curve. In such a model the asymmetry can become pos- 
itive at larger values of —t', caused by a negative real 
part in £ . The dash-dotted curve arises from an alter- 



native GPD approach 3J|, in which the imaginary part 
of 7i becomes negative while the real part of £ remains 
positive at larger values of — t'. 

An attempt to evaluate the complete set of Fourier 
amplitudes ||7J), and in particular the value of A^r^ - ^ , 
is presented in jl7j |. In this model, the GPDs are calcu- 
lated in a similar way as in the models 15j, [35j , except 
that the experimental value of the pion form factor F„ 
is used. Here a large non-pole contribution from E over- 
compensates the pion-pole contribution leading to the 
zero-crossing behavior of the amplitude as a function of 
— t' (see dashed curve in Fig. [2]). This model appears to 
be qualitatively in agreement with the data. However, 



(<p-<Ps 



IS 



within the large experimental uncertainty A UT 
also consistent with zero. A vanishing Fourier ampli- 
tude in this model implies the dominance (due to the 
pion pole) of E over H at low — t' . This is in agree- 
ment with the recent Hermes measurement of the ex- 
clusive 7r + cross section 221, which is well described at 



-t' = 0.1 GeV 2 by a GPD model |3| based only on E 
while neglecting the contribution of H. 

In summary, the Fourier amplitudes of the single-spin 
azimuthal asymmetry are measured in exclusive electro- 
production of 7r + mesons on transversely polarized pro- 
tons, for the first time. Within the experimental uncer- 
tainties the amplitude of the sin(</> — <pg) modulation is 
found to be consistent with zero, thus excluding a pure 
pion-pole contribution to the GPD E in leading-twist cal- 
culations. This could also be an indication for the dom- 
inance of E over the GPD H at low —t'. The observed 
amplitude of the sin 4>s modulation is large and posi- 
tive which implies the presence of a sizeable interference 
between contributions from longitudinal and transverse 
virtual photons. A next-to-leading twist calculation as 
well as knowledge of the contributions from transverse 
photons and their interference with longitudinal photons 
are required for a description of the measurements. 
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FIG. 2: Model predictions for the sm(<f> — 4>s) Fourier ampli- 
tude as a function of — i! . The curves represents predictions 
of GPD-model calculations. The full circles show the values 



of A 



sin(0-0 s ) 
UT,J 



taken from Fig. [T] The error bars (bands) rep- 
resent the statistical (systematic) uncertainties. See text for 
details. 
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